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Abstract. We have studied time variabihty in the flux from the flat spectrum syn- 
chrotron radiation of the Blandford & Konigl (1979) model for relativistic, conical jets. 
The resulting model has been applied to the flux variation of the flat spectrum of GRS 
1915+105 observed by Fender & Pooley (2000). This comparison has highlighted a funda- 
mental problem of the Blandford & Konigl (1979) model in that it requires unphysically 
large electron densities to explain the flux levels observed from the flat spectrum of GRS 
1915-K105. 



1. The Model 

To explain the infrared and millimetre emission oscillations from GRS 1915-1-105, 
we have adopted the Blandford & Konigl (1979, hereafter BK79)[0 model of 
a relativistic jet, containing power-law electron distributions. Adiabatic expan- 
sion gives the jet a conical shape with an electron density and magnetic field 
strength that decreases with radius. The jet emits partially self-absorbed syn- 
chrotron emission over a radial range covering several orders of magnitude, and 
hence from a wide range of optical depths. Therefore, the total jet emission con- 
sists of the summation self-absorbed synchrotron spectra from regions of the jet 
of decreasing optical depth (or increasing radius), and hence electrons at each 
radius, from the base to the top of the jet, emit a spectrum that is progressively 
shifted towards lower frequencies. This summation results in a 'flat' region to the 
spectrum, which has a slope of zero if energy losses due to adiabatic decompres- 
sion are ignored (as in the BK79 model). However, if some energy losses do occur 
then this region has a positive, inverted slope, the magnitude of which depends 
upon the magnitude of the energy loss. 

The extent of the flat spectral region, t'max/^'min) is solely determined by the 
emission region size, which in our model is fixed by the jet velocity and time- 
scale of variability. At the lower frequency end of the flat spectrum there is a 
smooth transition to an optically thick synchrotron spectrum with tx v^/"^, and 
the high frequency end smoothly terminates into an optically thin, F^, oc z^~^/®, 
synchrotron spectrum. 

Full detail of the derivation of our model may be found in Collins, Kaiser & 
Cox (2002)0]. The emission spectrum predicted by this model is given in terms 
of the optical depth function, Ty{r). 
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Dj is the distance to the jet, rmax/^'min defines the extent of the emission region, 
and all parameters have SI units unless stated otherwise. 

The model has two parameters. The first, which we have named ag, controls 
the optical depth as a function of radius, Ti,{r), and thus shifts the flat region 
of the spectrum to higher or lower frequencies. The second model parameter, 
which we have named 6o; only effects the flux normalisation. The ao parameter 
also has an effect upon the flux normalisation, and therefore we first fit the 
ao parameter to the spectral shape, before fitting 6o to the flux normalisation. 
The model parameter, oq, is related to the physical parameters of the system 
by ao = w{r = ro)/c(r = ro)B{r = ro)"^^/^, where k is the normalisation value of 
the electron density distribution, and B is the magnetic field strength. The bo 
parameter is defined as bo = w{r = ro)B(r = ro)^/^. The half- width of the jet 
w{r = ro) is fixed by observational constraints, and therefore from ao and 6o) we 
may determine k{r = ro) and B[r = ro). 

To implement time dependence in this model as a variation in the flux level 
of the flat spectrum it would be simplest to make 6o a function of time with 
ao fixed. However, physically this would require the electron density to vary by 
a process that exactly compensates variations in the magnetic field strength. 
Therefore we choose to make ao a function of time, and fix ho, which may be 
simply interpreted, physically, as a variation to the electron density injected into 
the jet. Time dependence is then included into equation |l| by defining rmax = v-^t, 
where v\ is the bulk velocity of the jet material, and by increasing the injected 
value oi ko with time according to a Gaussian function. 

From this model we can predict three observable parameters to which our 
model parameters may be fit; the time lag between the peak fluxes at each 
frequency of observation, the flux ratio of these peak fluxes, and the flux nor- 
malisation value. The value of the ao model parameter affects all three of these 
observables. The flux ratio is determined by the position of the flat spectral region 
with respect to the two frequencies of observation, which is defined by ao- The 
time lags are also determined by the position of the flat spectral region, because 
the optically thin end of the flat spectrum is formed at small radii, whereas the 
optically thick end is formed at much larger radii. Therefore the closer that both 
frequencies of observation are to the optically thin end, the smaller the time lag. 

2. Application to GRS 1915-|-105 

Fender & Pooley (2000) j^] observed large amplitude oscillations in the emission 
from GRS 1915-1-105 at an infrared wavelength (2.2 jim) and a millimetre wave- 
length (1.3 mm). The time lag between the peak fluxes from each wavelength for 
the first oscillation was approximately 25 s, with the dereddened infrared flux 
reaching 525 mJy and the millimetre flux reaching 340 mJy, giving a peak flux 
ratio of 1.5. However, the extinction correction applied to the infrared flux is 
uncertain, such that the peak flux ratio is 1.5 it 0.6. Such a small flux ratio over 
a frequency range covering three orders of magnitude is suggestive of the 'flat' 
spectrum produced by relativistic jets. 

For the fixed parameters of the model we assume that the distance to GRS 
1915-1-105 is 11 kpc, the opening half-angle of its jet is 4°, the base of the emission 
region is at rmin = 10^ m, and the jet material has a bulk velocity of = 0.6c. 
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Figure 1. Each line in this plot represents the 1.3 mm flux level versus time lag relationship for 
a different value of the peak injected electron density from the bottom line of fc(r = rmin) = 10^^ 
(solid line), 10^^, 10'^*, 10^*, to 10*° (dashed line) m~^ at the top. It is clear that for short time 
lags (from small values of ao) it is impossible to achieve the observed flux densities without 
extremely high values for the electron density. Taken from Collins et al. (2002). 

The value of the oq parameter was fit to the observed time lag, and, for this value 
of ao, ho was fit to the flux normalisation of the peak millimetre flux. From these 
model parameters the corresponding physical parameters of jet were determined 
to be S(r = rmin) = 7.8 x lO"'^ T and k{r = r^i^) = 5.4 x 10^° m'^. The model 
also predicts a flux ratio of 1.1 which agrees with the observations. 

We have found that it is impossible to obtain flux levels of the order of ~ 100 
mJy without unphysically high electron densities, for time lags of less than 300 s, 
as demonstrated by figure 1. This is not just a problem with this time dependent 
model, but is a fundamental problem of the BK79 model, as can be seen by 
performing the same fitting process to the observed flux ratio. 

To agree with the observed flux levels and the time lag between peaks with 
realistic electron densities, the model must either predict a flux ~ lO'* times larger 
with the fitted value of /cq , or predict the observed time lag with a value of oq that 
is ~ 10^ times larger. Since relativistic effects are negligible in the GRS 1915+105 
jet (see Collins et al. 2002), the later case can be immediately ruled out, as it 
would require the model's fixed parameters to be adjusted beyond physical limits. 
It is also difficult to justify the higher flux required by former solution, as Doppler 
beaming and geometrical effects are not sufficient. 

In conclusion it is impossible to explain synchrotron emission of this strength 
from a BK79 type jet with justifiable electron densities, when restricted to the 
observed time lags. 
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